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ABSTRACT 


An analytical approximation to the image theory model is developed for predicting the 
acoustical pressure field in a wedge-shaped ocean. This thesis is a continuation of the ongoing 
development of this model. Previously, limitations has restricted the source and receiver to be placed 
in only an upslope/downslope configuration. This thesis removes these limitations and allows the 


source and the receiver to be placed in cross-slope configurations. 
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I. INTRODUCTION 


Current historical events have necessitated a change from 
studying the propagation of sound in deep water to studying 
the propagation of sound in shallow water. With the fall of 
the Soviet Union and the rise of Mid-Eastern countries, 
Antisubmarine warfare (ASW) is more and more likely to occur 
in shallow water. Shallow water provides a much more complex 
environment for ASW prosecution, and, unfortunately, has not 
had as much research as in deep water. 

Image theory is one method that is being developed to 
predict sound propagation in shallow water. This method is 
more accurate than ray tracing methods because it is not 
limited to high frequencies. A number of thesis conducted at 
the Naval Post Graduate School were involved with applying 
image eneery to a wedge-shaped, shallow-water ocean. The 
first thesis provided the groundwork for the initial 
development of the image theory. The most recent 
[Nassopoulos, 1990} derives a simplification involving a 
doublet approximation and incorporates all of this information 
into one manageable computer program [Nassopoulos, 1992}. 

This thesis assimilates a more complicated cross~slope 


approximation into the doublet analysis which allows the 


source and receiver to be placed anywhere in the wedge. A 





computer program has also been developed to show that the 


approximations are correct within the set of assumptions. 





II. BACKGROUND 


Research on the prediction of sound propagation in a 
wedge-shaped ocean has been directed in three major areas: the 
parabolic equation approximation, the adiabatic normal mode 
theory, and the image theory. 

The parabolic equation can be used as a range dependent 
underwater acoustic propagation model. It basically replaces 
the Helmholz equation with a one-way parabolic approximation 
which will generate the acoustic field as an initial value 
problem. The original version of the parabolic equation, 
developed in 1977 by Tappert, contained restrictions on the 
maximum elevation angle at the source. These limitations were 
reduced by Claerbout and Green and finally eliminated with the 
use of a higher order parabolic equation (Collins, 1987]. 
Because the parabolic equation is cylindrically symmetric, 
sloping bottoms cause a problem. One attempt to solve this 
problem has been attempted by sectioning off the bottom into 
a series of range-independent regions and then applying the 
normal parabolic approximation to each section {Lee et al, 
1983]. Although this method increases the accuracy, it still 
does not alleviate the problem completely. A more recent 


approach to this problem has been attempted by using a rotated 


parabolic equation [Collins, 1990] which shows good results. 





Normal mode theory is a range-independent approach. This 
approach has been used for open ocean predictions, but for a 
wedge shaped ocean, a range dependent approach is needed. For 
this type of approach it was shown [Pierce, 1965] that an 
approximation to the normal mode theory can be utilized by 
performing an adiabatic separation of the depth and range 
coordinates in the wave equation. This approach has been 
applied to the wedge-shaped ocean [Graves et al, 1975] but was 
limited to only small slopes and non-penetrable bottoms. An 
exact solution for the ideal wedge, where boundaries are 
perfectly reflected, was developed. [Buckingham, 1987] A 
continuing problem for adiabatic mode theory is that it can 
not explain the transition from the guided mode to the 
evanescent modes at cutoff (Jensen et al, 1980]. One hybrid 
approach that has been pursued to account for this is the use 
of the parabolic equation [Pierce, 1982] which is solved 
asymptotically subject to adiabatic mode initial conditions. 
Another method is based on tracking of local plane-wave 
spectra [Arnold et al, 1984] which produce intrinsic mode 
fields. Another method uses spectral representation based 
on the theory of characteristic Green’s functions {Kamel et 
al, 1983]. 

Another approach in solving the wedge-shaped ocean problem 
is image theory. This theory predicts the pressure field from 


a series of projected images [Coppens et al, 1980]. A further 


description of this approach will be provided in the following 








rw 


chapters. The basic geometry and mathematical analysis has 
been developed and compared to the parabolic equation model 
with favorable results. The ups lope/downslope case has been 
approximated [Nassopoulos, 1992] by use of Taylor Series 
approximations resulting in a substantial savings in computer 
computational time. Although less research time has been 
placed into the image theory approach, it shows high potential 
as a viable way to predict acoustic pressure fields in a wedge 


shaped ocean. 


III. DEVELOPMENT 


The basic goal of this research is to approximate the 
images by a sum of doublets. Although a full description of 
the theories will not be incorporated into this thesis, the 


following can be used as a general overview. 


As. IMAGE THEORY 

Image theory stems from the idea that if a source is 
radiating near a surface, then an image, out of phase with the 
source and an equal distance away from the surface, is 
produced. In an isospeed wedge-shaped ocean, numerous images 
are produced [MacPherson et al, 1966}. These images form a 
circular pattern around the point where the surface meets the 
bottom and interact with a series of images each possessing an 
image of the surface and bottom of the wedge. The images are 
numbered sequentially around the upper half of the circle 
starting with 1 for the source, 2 for the image and so on 
around the circle. The lower half of the circle is labeled in 


a similar fashion. The pattern for the sign convention along 


with the numbering of the images is shown in Fig. 1. 





Surface 


Bottom 





Image Pattern for a Wedge-Shaped Ocean 
Figure 1 
Each image contributes to the total pressure field at a point 
in the wedge. For each image, the distance to the receiver 
and the total reflection coefficient, which is the product of 
all the reflection coefficients along the path from the image 
to the receiver, can be determined. The total pressure is the 
sum of the pressures from all of the images. A simplified 


three dimensional view at Fig. 1 shows that the shoreline acts 


as the perpendicular axis to the circle of images. See Fig. 





3-Dimensional View of Wedge 
Figure 2 
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B. DOUBLET RADIATION 

The images on each side of the pressure release surface 
can be grouped together in pairs to form acoustic doublets: 
Images 1 and 2 form a doublet, 3 and 4 form another doublet, 
and so on around the upper and lower half of the circle. Each 
soubiee consists of two simple sources, vibrating at the same 


frequency but 180 degrees out of phase. See Fig. 3. 





P (8,7, t) 


Geometrical Presentation of Doublet 
Figure 3 


It can be shown that in the far field the pressure field 


resulting from this doublet using a far field approximation is 


P(r, 0, t)=jAsin( 2 kdsin®) ef (ween (a) 


Where ad is the separation of the source and image in the 


doublet, r is the distance from the field point to the 


midpoint of the doublet, and @ is the angle formed by r and a 
perpendicular bisector of the doublet. It should be noted 
that, if the polarities were switched, the pressure field 
would change sign. 

To incorporate doublets into the wedge problem, Eq. 1 
needs to be modified to include the reflection coefficients. 
Nassapolous developed the following equation for pressure 

A 1 
P(r,6,t)=—=[(x-) cos (=k, dsin@) 
I ; 2 (2) 
+3 (x+y) sin(>k,dsin6) ] e7 (we-kr) 
Where x and y are the cumulative reflection coefficients for 


the negative and positive images respectively. 


C. DOUBLET ANALYSIS OF THE IMAGE THEORY 
1. General Description: 

The purpose of this thesis is to extend Nassapolous’ 
doublet approximation to the image theory model to ieakoonace 
cross~slope propagation. The basic assumptions remain the 
sane the sound speeds and densities are constant in both the 
wedge and the bottom, the water-air boundary is a pressure- 
release surface, and the slope is constant and shallow 
(Nassapolous, 1992]. Nassapolous developed the case for 
upslope and downslope propagation. Figure 4 shows how 
doublets are grouped together: The third and fourth upper 


images are grouped together with the first and second lower 


images (n=1), the fifth and sixth upper images are grouped 





together with the third and fourth lower images (n=2) and so 


on around the circle until all the doublets are paired. 





Representation of the Groupings of Acoustical Doublets 
Figure 4 


The source and image for the first surface reflected path form 
the neutral doublet (n=0) and is an exception to this grouping 
rule. Similar to doublet radiation, where each doublet 
contributes to the pressure field, now, each pair of doublets 
is considered the basic contributor to the total pressure 
field. The number of paired doublets used in this total 


pressure calculation can be determined by 


(3) 


B is the angle between the surface and the bottom. To ensure 
that the pressure field is accurate, a wedge angle should be 
chosen such that all of the images are grouped into pairs of 


doublets. To insure that there is closure for the set of 
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images, 360/@8 must be a multiple of four. In calculating the 
total pressure field, the pressure from the neutral doublet 
‘will first be developed, followed by the pressure field from 
each of the paired doublets. 
All distances used are scaled according to Ryag 
=R,.3/X- Two bottom types are examined. in this research, a 
fast bottom and a slow bottom. For a fast bottom, the scaling 
distance X, occurs where the lowest mode attains cut off. For 
a slow bottom, the lowest mode does not attain cut off, Hue a 
scaling distance X, is used for convenience in calculations. 
The definitions of these scaling distances are 


Fast bottom: 


x Eiper eh ct alt 
ws 28inO.tanb 
T 
Ke xX 2a ee 
a"* 2tan@ tanp (4) 
C, 


cos8,=— 


N 


Siow bottom: ‘i 


kK, X= ———_ > 
ws 2tané,tanp 
T 


KX = 
"5 2s5in6,tanB +3) 


C. 
cos8,=— 
a 
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B is the wedge angle, the subscript 1 refers to the fluid in 
the water, and the subscript 2 refers to the bottom 
{Nassopoulos, 1992]. 
2. Neutral Doublet (n=0) 
The neutral doublet, as stated above, comprises of the 
source and upper image. The pressure field from the neutral 
doublet crosses no bottom boundaries therefore, no reflection 


coefficients need to be considered. See Fig. 5. 





Geomerrical Representation of Neutral Doubler 
Figure 5 


d is the distance between image 1 (the source) and image 2, 6 

is the angle formed between the surface and the plane 
containing the receiver and the shoreline, and R, is the 
perpendicular distance from the shoreline to the receiver. 
‘Taking Eq. 1, that was developed for doublet radiation without 
reflection coefficients, and applying it to the new cross 


slope geometry gives 


PJ-j 22 sin ( $k,dsino,)) 


oO 


ad (wt-kRo) 


(6) 
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where P’, is the pressure field from the neutral doublet, R’, 
is the distance between the receiver and the midpoint between 
the neutral doublet and oa’, is the angle formed between R’, and 


the surface. 


To reduce Eq. 6, the following substitution is made: 


$-R,siny (7) 


y is the angle formed between the surface plane and the plane 


containing the source and shoreline, and R, is’ the 


perpendicular distance from the shoreline to the source. See 


Fig. 6. a os 


Surface 


Source 


Geometrical Representation of R; and ¥ 
Figure 6 : 


Since 6 is assumed to be a small angle, siny can be 


approximated by y which reduces Eq. 6 to 


; . a7 eae a 
pJ=j24sin(k,Rysino,)) ef (8) 


° 


A further substitution can be made by applying the geometry in 


Fig. 7. 
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Receiver 


Gcornettical Representation of —) and R’y 
Figure 7 


The side common to both triangles shows that R,siné = 


R’,sino’ Since 6 is small, the small angle approximation can 


be utilized to form the equation sino’, = R,6/R’,. 
Substituting this into Eq. 8 and assigning [’, = k,R,R,/R’, 


gives the final pressure for the neutral doublet 


p,/=j #4sin(T/8) e 
Re 


j (we-k,R,’) 


(9) 


3. Pressure Field From the Nth Pair of Acoustic Doublets 

The paired doublets can be split into two groups, one 

from the upper set of doublets P’,* and one from the lower set 
of doublets P’,. The "'" denotes the cross slope geometry, 
the subscript n denotes which paired doublet is being 
calculated and the "+" and the "-" denotes the upper and lower 
doublet respectively. Figure 8 illustrates the geometry of 


these upper and lower doublets. 
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U 
doublet 





nth 
lower 
doublet ; 





Geometry of Upper and Lower Doublets 
Figure 8 
The pressure from the upper doublet is 





Pft=(-1) 724 sin(ik,dsing,!*) ei wee) (10) 
Re 2 
and the pressure from the lower doublet is 
. oF, f- 
Pp! =-(-1)°7 24 sin(+k,dsino,/”) ene (11) : 
R,, 
The (-1) and the -(-1) account for the different 


orientations of the top and bottom doublets. 


The total pressure from all of the paired doublets is 


N 
P/=PJ+)) PLtP, (12) 
n=1 


-This pressure could be calculated by a computer program, but 
it would take a relatively long time to compute. To shorten 
the computational time, and also research for an analytical 
approximation a Taylor series expansion is used for R’,’, 


R’,, Sino’,*, and sino’,. 
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a. Derivation R’,” and R’,- 
The geometry for an upper doublet image is shown 


in Fig. 9. 







Receiver 
Plane 


Receiver 
aa ey ee] 
Geometrical Representation of Upper Doublet 
From this geometry Figured, 
(RI) teh) es (y)* . (13) 
and 
(R,,*)?=R,2+R,?-2R,R,cos (2nBp +8) (14) 


therefore by direct substitution 


Ri*= R,°+R,*+Y,?-2R,R,cos (2nB +6) (15) 


By the same reasoning, for the lower doublet 





R, =/R,?+R,?+Y,7-2R|R,cos (anp-6) se 


Setting Y, to zero, Eq. 15 and Eq. 16 reduce to 


R, =R,2+R,?-2R,R,COS (2np +0). (17) 


~=/R,?+R.7-2R,R,cos (2np-8) (18) 


which match the equations developed be Nassapolous. 
(1) The Taylor Series Expansion for R’,” and R’,. 


Because R’,” and R’, appears in the phase a third order Taylor 
series expansion is needed. Taking this expansion around 6=0 
gives 

3p/ 
2 aR, 1 93 A" Rn 


“ype 80% 37 Oye |o-0 oa 





1 
Ri=Ri +b Oe = 2 hot 3° 


This formula is the general case for both the upper and lower 


doublet. To get the correct formula, R’, would be replaced 


n 


with R’,” or R’, as needed. R’,, is the distance to the nth 


no 


doublet if the receiver were moved vertically to the surface. 


See Fig. 10. Since the geometry is now symmetrical 


Rio =Rio =Rjo=fRy2 +R, + ¥,°-2R,R, COs (2p) (20) 








and Rng 


Geometrical Representation of ¢ oe 
Figure 10 


Taking the first derivative yields 


dR, _ R,R,sin(2np+5) 
dé RO 





dR, -R,R,sin(2np-8) 








dé Re 
dri ae R,R,sin (2nB) 
db Ro 
dR}, — ~R,R,sin(2nB) 
dé '° R! 
no 


Taking the second derivative yields 
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(22) 


(22) 


(23) 


(24) 








d?Ry _ [Rn R,R,cos (2nB +8) ] 





dé? me 
a (25) 
_ [R,R,8in(2nB +6) R,R,sin (2nP+5) J 
(z.)? 
d?R,” R,R,cos(2nB)  R,?R,?sin? (2n6) 
—™ |,..25 2 eee (26) 
db? Fino Re 
2plr RR +2 
IRS ge AF [cos (2np) - asin {2aP)) aa) 
dé Ro Rno 


Through cancellations of the minus signs, the second 


derivative for the lower doublet is the same, 


so 
d@Ri* d?*R, (28) 
ba 8 
dé dé 
Taking the third derivative and evaluating at 6=0 yields 
dR," R,R, , ~3&,R,cos (2nB) sin(2n | 
ied ae ce 1%, (2nB) (2nB) _ sin (2np) 
dé} pe / 2 
no Rno (29) 
. 3R,?R,' sin? (2nB) , 
74 
Rno 
d?R, R,R -3R,R,cos(2n in(2n 
as | 9 lect y a8 a COs 2 SE) tee oApl 
ab° Ri Re 
no (30) 
P 3R,7R,7sin* (2np) ; 
4 
Rao 


Combining these pieces together, the third order Taylor series 


expansion is obtained 
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R,R,sin(2np) R.R 
/ 1 
rag = Riot a + 58 [cos (2nB) 


Rno no 
7 caeeee enh) kek RR 83 | ~3R,R,cos (2nB) sin(2nB) (32) 
Rno no RS - 
. SRR sin? 
- sin(2nB) + Saree RAE 

14 

Eno 
(2) Phase Approximation. R’,, and R’, play 
significant roles in the phase. Placing these terms in the 


phase in the third order Taylor series expansion and 


rearranging the terms gives 


R,R,sin (2nB) , BPs 


cxp (Jwe) exp $ (JB) [ (Rt) + 6 ——- 6? cos (2nfB) 


7 
Rno 2Rno 
_ R,R,sin? (2nB) ~ Fi®ags -3R,R,cos (2nB) sin(2nB) 
Ree 6 Rno Ro 


2p 2a47n3 
-nintaapy + 2eteteP eet) T) 


no 
(32) 
Since e™ does not contain the index n, it can be set aside for 
the moment and the second term will be examined for possible 
reductions. In order for a term to be considered negligible, 
a phase error of less than 7/4 is assumed. Certain 
assumptions also need to be made in order to represent a 


typical ocean environment. First, it is assumed that 


J |<0.1 (33) 
1 


which is equivalent to 





sin0.<0.4 


(34) 


Where c, is the speed of sound in the water, c, is the speed of 


sound in the bottom, and @, is the critical angle. 





Second, 
B<0.1 
Finally, 
Sigel 
Re 10 


(35) 


(36) 


Whose implementation will be seen in the following analysis. 


Grouping the terms inside the brackets from Eq. 32 in 


terms of powers in R, and R, gives 
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an [Ssin’ (2nB) 83} (37) 

PR 1 gin? (2np) 8*-2.cos (2np) sin (2nB) 8") (38) 
13 72 2 

k-2e (sin (2np) b+ 3.008 (2nB) 6 sin (208) 8° (39) 

kKRi, (40) 


Each group needs to be compared to 7/4 to determine if it can 
be considered negligible. Beginning with expression 37 the 


following inequality is set up 








Ry Red oe wes 
k “2 [>sin* (2nB) 6 je«n/4 (41) 
nO 
Rearranging gives 
R, 2 Rk, @ 4. + 43 3 
k=) (—-) R(>sin (2nB) 6?) <«n/4 _ (42) 
NO NO 


Making the definitions S=R,/R’,, and sinn=R,/R’,, and multiplying 


through by X/X gives 


kX s°sin?Q—2 (sin? (anp) 8°]<n/4 (43) 


where kX=7/(2sin@,tan8) for a fast bottom. 


To ensure that the left hand side of the equation is always 








<< 1/4 the maximum of all powers of sine and cosine terms will 
be taken to be unity, with the exception of sin§#, which, as 
stated in the above assumptions, already has a maximum value 
of 0.4. Substituting these maximum values into Eq. 43 and 


rearranging gives 


gs ke (44) 
.4tanB” XxX ~ 83 


Since 8 is small, tan 6 = 6B while &, the receiver angle, falls 
between 0 and 6. The right hand side is minimized by setting 
B=6. Making this substitution into Eq. 44 and allowing B to 


take on a minimum value of 0.1, reduces Eq. 44 to 


R 
g32<40 (45) 
x 


Finally, substituting S = 1/10 according to Eq. 36 gives 


R 
—2<40000 (46) 
x 


If R,/X for some reason became larger than 40000, the term 
would need to be added back into the equation to keep the 
phase error less than 7/4. Note that 6 2 0.1 is eqivalent to 


setting a lower limit on the wedge angle of about 6 degrees. 


Expression 38 can be written 





R R 
k(—)?( a )R, [Ssin? (2nB) 8 


2 
no no 








(47) 
- cos (2nB) sin (2nB) 8°] «n/4 


Making the same substitutions as in the previous equation 


yields 


t 2 Re +2 2 
psibeenh (48) 


- cos (2nB) sin (2nB) 67] «x/2 


Once again substituting in the maximizing values gives 


g2 ft eg (49) 
XxX 


Substituting S = 1/10 gives 

R : 

—2<400 (50) 
x 


which allows this expression to be neglected. 


Expression 39 gives 





/ 
nO 


“1p, [sin (2nB) 8 +cos (2n6) §¢-Zsin(2nB) b°1<7 (51) 


Making the same substitutions as in the previous equation 


yields 
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Te 


Ree. 4 1 ; 
Zambccan, et eth coat ane = 


1 (52) 
| sin (2nB) 8?) <i 


Rearranging the terms and substituting in the maximum values 


gives 


R 
S—e« 





02 
1 


53 
sin(2nB)& + Scos (2nf) 8? ~ =sin(2nB) d° oe 


The dominate terms in the denominator are sin(2n8)6 + 
0.5cos(2n8)6’. Since 1/6 sin(2n@) 63 is a comparatively small 
number it can be neglected without increasing the phase error 
over 7/4. 


It follows that the expression for the upper doublet is 








“Jk (Rgo* “272 (sin (2np) 8+ bcos (2nf) 6) } (54) 
e no 
and by the same method 
~jk(R/,- “2% (gin (2np) +2cos (2nf) 8?) } 
no R! 2 (55) 
e no 
| for the lower doublet. 
| 
3) Amplitude Approximation. It can be shown 
1g) PP 


that, in the far field, the amplitude, R’,* and R’, can be 


replaced by R’,, with a relatively small margin of error. 
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Gj ¢ 2 7 + ‘ om 
b. Derivation of sino’,’ and sing 


n 


and 


a 


To take the Taylor series expansion of sino’,* 


sino’,, an expression for each term is needed. The geometry 


nf 


for this is shown in Fig. 11 for an upper doublet. 


Apex 





Receiver i 


Geometrical Representation of g++ andR,* 
Figure 11 ; 


From the law of sines 


ae 
Ra _ VR * Yo" (56) 


i ‘ / 
sinn sino), 





Rearranging this gives 
Ri’ sino) =sinn/R,?+Y, (57) 
Substitution WV1-cos’ 7 for sin 9 gives 


R, sing), =/(R,2+¥,7) - (R,2+¥,7) cos? (58) 


To get » in terms of 6 a substitution for cos’ 7 is needed. 


To do this the following equality is setup: 





+2 3 
Ry =R,?+R,?+Y,2-2R,R,cos (2nB+8) = 


(59) 
Ri +R7+Y5-2R,fRF+¥2cosn 
Canceling similar terms and solving for cos? y gives 
Récos? (2n 
pogine neces ea) 70): (60) 


R2+¥?2 


Substituting this into Eq. 58 and replacing cos*(2n8 + 6) 


with i-sin’(2n@ + 6) gives: 


Ri’ sinoy’=/Y,7+R,’sin? (2np+d) (62) 


Bringing R’,” over to the right side gives the final 


+ 


expression for sino’,*. 


one f¥ 2 +R, Sin? +6 
sino, _V¥o° +R,*sin* (2nB +6) (62) 


= = 


Rp 


And finally, by the same reasoning the expression for sing’, 


is 


, y-_ Vf ¥,°+R,*sin? (2np-6 
sino), ea ee (63) 
Rn 


When the receiver is moved vertically to the surface Eq. 62 


and Eq. 63 become 


; VY,°+R,’sin* (2n 
sino’ Vo Ha San UP) (64) 


no” / 


Rao 


Again, by setting Y,=0 for a geometrical check 
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+_ &,Sin(2nB +8) 


sing 
n R,* (65) 
. _ R,sin(2np-8 
sing, = Zap) (66) 


which agrees with the upslope/downslope case. 


(1) The Taylor series expansion for sino’,” and 
Sino’, . The first order Taylor series expansion taken around 


an 


Y,=0 gives 


sine,/=sing,./+8 0 |,,, (67) 
Because these terms only appear in the amplitude, the 
approximation sino’, * sino’,, is acceptably accurate. 
c. Total Pressure Including Approximations 
Putting all of these equations together yields the 


pressure for the upper doublet, 








Pj’ =(-1) 2724 [sin>k,dsinon,] ei (wt) 
Pno (68) 
-jk, [RAg+ 23 (sin (2nB) 6+ 5 cos (2nB) 8?) | 
e nO 


That for the lower doublet is the same except for a sign, 
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Pi=s-(-1)3 eee A tsini has .dsino!, Jed(we) 


0 (69) 
’ /_ RR, : 1 
~jk, {Rao 7 (sin(2nB)$+=cos (2nB) 8?) j 


e ‘2O 





The pressure from the pair is 


ee 1) 2 SA 4 [sin= 1k dsino), 5] ed (WEE Rao) 
rR 
no (70) 
k,R,R : 
sin{—++# (sin (2nB) 8+Scos (2nB) 82) ] 


no 


The total pressure as stated in Eq. 12 is 


N 
P/=P J+) PtP, (71) 
n=1 


Combining Eq. 9 and Eq. 70 gives the final solution for the 


total pressure, 


P,=j #4 sin (Toy) e 
R! 


/ 
wae 1)" [sin Sk _dsino!,] e7 Ye Fa! (72) 


no 


Pal wt-k,R4) 3 


sin [| 2% (sin (2nB)6+4cos (anf) 8°) 
Rno 
4. Validation 
To validate the analysis of doublet cross-slope 
propagation, a computer program was developed. See Appendix 
"A", Different inputs of Y,, R,, and R, were used along with 
a wedge angle and source depth of 7.5 and 3.75 degrees 
respectively. These inputs were also put into the program 


URTEXT anda comparative table was constructed. See Table 1. 
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TABLE 1 


COMPARISON OF DOUBLET ANALYSIS TO IMAGE THEORY 
p=7.5 y=3.75 # Of depth points = 6 






Ri=1 R2=15 YO=0 Rl=1 R2=12 YO=8 
DOUBLET _. IMAGE DOUBLET IMAGE 


| AMP | PHS | AMP | PHS [ AMP | pHs | AMP | PHS 
99135 
99132 
40713 | 99133 
49862 | 99135 
3 | -1.203 1.180 99135 
99138 
AAMP = .0066_APHS = .01898 

Rl=1 R2=14 YO=3 [| Ri=1 R2=8 YO=9 | 

DOUBLET IMAGE 
PHS 


7774 


.31909 | .7759 





ann 

a 

> 

Oo 

a | 
Olu 
Go| 
an 
mil 
alt 
O1O 
eat ie-) 
0100 
w]e 
cnt oo 
WG 
IO 
AE 
.! .! 
NIN 
MO 
IN 
onl Kon 
Nil 
uty 
an) Ie) 
Ol h. 
Ci} 
.! ! 
NIN 
bib 
al fo 
NO} Go 










4 
5 
42229 | 80677 
51721 | 80558 
AAMP = .0423 APHS = .0351 


Rl=1 R2=13 YO=6 1=1 R2=2 YO=10 


DOUBLET IMAGE DOUBLET IM 


w 


> 
oO 
tT 





76982 | .40469 
76894 
71336 






AAMP = .11073 APHS = .0272 AAMP = .0031  APHS = .3421 











It is shown that the doublet cross-slope propagation results 
are in good agreement with the URTEXT results. It should be 
noted though, that for the results to compare favorably, the 
wedge angle that is chosen must ensure that the number of 


wedges are an even multiple of four. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


Image theory is an excellent method for predicting the 
transmission of sound in shallow water. It is an important 
area of research because it is directly applicable to shallow 
water Anti-submarine Warfare. Other methods that are being 
examined include the Parabolic Equation approximation, ray 
tracing methods, and the Adiabatic Normal Mode theory. Image 
theory shows the most promise for predicting the acoustical 
pressure fields in a wedge-shaped ocean because it can predict 
pressure fields for cross-slope and for all frequencies. 
Image theory can also deal with the transitional area at cut 
see ee region that is difficult for normal mode theory. 

Previous theses have laid the groundwork for the image 
theory model. This thesis has extended the doublet 
approximation to cross-slope propagation. The base 
approximation used, R,? << R,’ + Y,’, limits the extension to the 
far-field. Another limiting restriction on this thesis 
basically states that for a phase error of 7/4 or less, the 
scaled receiver distance» R,/X must be less than 400. 
Geometrically this says that when the receiver gets too far 


away the phase errors combine and the pressure equation 


developed in Eq. 72 is no longer valid. One further 





restriction is that the. wedge angle must be chosen such that 
the number of wedges produced is an even multiple of four. 
Because of this, when the dipoles are grouped into pairs of 
doublets, one dipole, 180 degrees from the neutral doublet, 
will be left ungrouped. This will only be a significant 
factor when the grouped dipoles cancel completely leaving only 
the dipole to contribute to the pressure. The effects of this 
dipole can also be seen in Table 1: The results vary when R, 
gets small and Y, gets large. Although this dipole can have 
noticeable effects on the pressure in the rigid bottom case, 
when reflection coefficients are incorporated, the effects 
become negligible. 

The next step in this research would be to incorporate the 
reflection coefficients for the upper and lower doublet into 
Eq. 72. The reflection coefficients though, because of their 
large effect on the pressure, need to be calculated 
individually for each image. Nassopoulos incorporated the 
reflection coefficients into the upslope/downslope case, but 


an extension is still needed for the cross-slope case. 
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APPENDIX "A" 


REM WHOSE DERRIE. DIROLE APPROXIMATION 

PEM crons slome. wedge with rigid bottom, including third order approx 
PTS 14 oOeH 

INPUT “analee( dec) of wedqe. source: B,G =",B,G 

INPUT "M1 = B/C delta O) =".M1 

NETS =INICCC INTE 360 “B+ 9001 )-2 )74 + 0001 ) 

INPUT “source, receiver. apex distances: R1L,.R2,YO =",R1.R2,YO 
INPUT “et7e2 —- "Ve 
T6=1B80/PI: B=B/T6: G=G/T4S 

STARCOS = SORCARSCC*C-1))! T4A=PI/( 2xSTARCOS*TAN( B)) 
LPRINT USING “THESIS APPROX angles: wedge= ##.4# source= ##.#°5B*T6 .G*T6 
LPRINT USING “distances: source=t#.## receiver=#t.## shore=kt.#"FR1.R2.YO 
LERINT USING “ cl/o2= #.#H "3C 

LPRPINT WING "kive=pi/(2 tanR sqrabs(1-c72) )=Hat  #eHe "S74 
LERINT USING "No. of Dipole-Pairs =#8.##";SNPTS 

LFRING * DEG ANP PHASE " 

DIST=Z2¢R1I¥SING G) 

YZRYORYO? R22eR2*H2 

S2=C*C2 DZ=VOXYORRIFRLIRZ*R2Z 

MR=2kREERZ 

FOR M=1 TO LOO: DO tR/M1i: IF (B-D)<-.0001 THEN 380 
ROZSARC D2-RB*COS(D))}! GAMOFRL¥*¥R2/RO 

Te( 2/ROIVSING TAXGAMO*GYTD ) 

PLET*SIN( T4*RO ): P2=T*COS( T4*RO ) 

P3asO: FARO 

FOR N#1 10 NPTS ‘ 
PN=SQR( P2-R3A*COS( 2 *NYB )) 

RNL=( T4*R1L*R2*SIN( 24*N*B ) )/RN 

RN2=( T4AkRL¥R2*COS( 2*N FR) )/RN 

SNS=( SARC Y2+R22*( SINC 2*N*B ) 72) )/RN 

Xa( Cet YN KC 2/RN SING .S*TATDIST#SNS ) 

THOPOS= ~TA¥RN - D*RNL ~ .5*072*RN2 

THDNEG= -T4*RN + D*RN1 -  .5*D72*RN2 

THOR = COS( THDPOS )-COSC THDNEG ) 

THD = SINC THDPOS )-SINC THDNEG } 

REALL = -X*THDI 

IMAG1 = X*THDR 

PAsREALL + 2 

PATIMAGL & 4 

NEXT N 

PREPLtP a: PIMPQEPA 

AMP=SOAR( PRIPRIPIM*FIM ) 

PHS=ATNOCFIM/PR ) 

LPRINT USING ° HHH HLTH 6 #4 HAHAH OS OKTE LAMP PHS 
NEXT M 

END 
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